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Using a semiclassical mean field theory, we show that the screening potential exhibits a 
characteristic radial variation in the tunneling region in sharp contrast to the assumption 
of the constant shift in all previous works. Also, we show that the explicit treatment 
of the tunneling region gives a larger screening energy than that in the conventional 
approach, which studies the time evolution only in the classical region and estimates 
the screening energy from the screening potential at the external classical turning point. 

This modification becomes important if the electronic state is not a single adiabatic state 
at the external turning point either by pre-tunneling transitions of the electronic state 
or by the symmetry of the system even if there is no essential change with the electronic 
state in the tunneling region. 

Nuclear reaction rates at low energies play the key role in energy generation in stars and the primordial 
and stellar nucleosynthesis. The bare reaction rates are modified in stars by the screening effects of free 
and bound electrons. The knowledge of the bare nuclear reaction rates at low energies is important not 
only for the understanding of various astrophysical nuclear problems, but also for assessing the effects 
of host material in low energy nuclear fusion reactions in matter. This is currently a subject of great 
interest in nuclear physics. Rolfs and his colleagues have reported that the experimental cross sections of 
the ^He(d,p)"'^He and of D(^He,p) ^He reactions with gas target show an increasing enhancement with 
decreasing bombarding energy with respect to the values obtained by extrapolating from the data at 
high energies [1]. They also claimed that the enhancement is larger in the ^He(d,p)^He reaction. Since 
then similar enhancement has been reported for many systems with not only gas targets, but also with 
metal targets such as the ®Li(p,a)^He reaction. 

These observations have motivated many theoretical as well as experimental studies. Many of them 
attempted to attribute the enhancement of the reaction rate to the screening effects by bound target 
electrons. A simple approach is to assume that the screening effects can be well represented by a constant, 
i.e. radially independent, decrease of the barrier height in the tunneling region. This decrease is named 
the screening energy. It is determined by making a fit to the data. A puzzle is that the screening 
energy obtained by this procedure exceeds the value in the so called adiabatic limit, which is given by 
the difference of the binding energies in the united atom and in the target atom and is theoretically 
thought to provide the maximum screening energy, for all systems so far studied experimentally [2] (see 
ref. [3] for a recent modification). For ^Li(p,Q!)a reaction, in addition to the direct measurement, an 
indirect measurement of the cross section using the Trojan horse method has recently been made [4]. 
The comparison between the two methods indicates again that the screening energy in the direct method 
exceeds the adiabatic limit by a large factor. 

One should, however, keep in mind that the stopping power is not well established, especially for gas 
target, at such low energies [5-7,3]. Also different values of the screening energy are obtained depending 
on the method of analysis [8,9]. 
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In this paper, we discuss the properties of the screening potential in the tunneling region. We examine, 
in particular, whether it can be represented by a constant shift as has been postulated in all previous 
studies. We also examine the validity of the former dynamical approach in refs. [10,11], which solves 
the coupled equations for the electronic and nuclear motions only in the classical region, and estimates 
the screening energy by using the electronic wave function at the external classical turning point. To 
that end, we describe the time evolution of the electrons by a Schrodinger equation and the relative 
motion between the projectile and target nuclei by classical Newtonian equations. They are coupled to 
each other through a variational principle leading to a mean field theory. In that sense, our formalism 
is the same as that of Shoppa et al. [10,11] for the classical region. However, we extend the study to the 
tunneling region as well. 

We denote the coordinate of the relative motion between the projectile and target nuclei by R and 
that of the electrons by which contains in general the coordinate of the center of mass of electrons 
relative to the center of mass of the target and projectile nuclei, as well as their intrinsic coordinates. 
Considering the head on collision, we assume the following Hamiltonian for the total system, 

where (i?) is the bare interaction between the target and projectile nuclei, Hq is the unperturbed 
Hamiltonian of the electrons, and Vc{R,^) is the interaction between the electrons and nuclei. Denoting 
the wave function of electrons and the distance between the projectile and the target at time t by t) 
and respectively, the time dependent Schrodinger equation and the classical Newtonian equation 

for them read 

= [^o(0 + 0] t) (2) 

{R) + AH(i?)] (3) 

where 

AH(R) = {4>\ 6] 1^) (4) 

Eqs.(2) and (3) lead to the following energy conservation law. 

Yi^f + y^°KRit)) + ^viR) = E ( 5 ) 

Eqs. (2) through (4) have been derived from the time dependent Schrodinger equation for the total 
wave function by approximating '^{R,^,t) by a product of the wave functions for the relative 

motion between nuclei t) and that for electrons and determining them by a variational 

principle, thus the name of a mean field theory, and then converting the Schrodinger equation for xiR^ R) 
into Newtonian equation for R{t) following the idea of Ehrenfest, i.e. by assuming that the classical 
time dependent variable R{t) and its conjugate momentum P{t) are given by the expectation values of 
the corresponding operators in the state x(t) (see ref. [12] for details). We call AV{R) the screening 
potential, which is nothing but the total energy of the electrons at each time t or equivalently at each 
R{t). We note that eq.(3) is equivalent to 

which is more familiar in literatures. We prefer to the expression in eq.(3), because it connects more 
directly to the corresponding equation in the tunneling region which we show below. 

We determine the time evolution in the classically allowed region by solving eqs.(2) and (3) along the 
real time axis with the proper initial condition. Once the velocity of the relative motion becomes zero, 
we switch to the imaginary time, t = —it, and continue to follow the time evolution in the tunneling 
region using the following equations. 
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[Ho{0 + Vc{R{r),0]<P{^,r) 


(7) 


M 


^ dr 

d^Rir) 


d 


V‘^^\R) + AV{R) 


dr^ dR L’ 

The screening potential and the energy conservation law in the tunneling region are given by, 

[i7o(?) + V;(i?(r),0]|</>) 


AV{R) = 


{<l^\4>) 


( 8 ) 

(9) 


+ v'^^\r(t)) + AV(R)=E ( 10 ) 

2 OT 

We note that the norm of the wave function of electrons is not conserved in the tunneling region. 
Accordingly, the denominator of the screening potential given by the r.h.s. of eq.(9) is essential as we 
see later. We note also that the potential renormalization given by eq.(9) is the equivalent potential of 
the dynamical norm factor, which has been introduced in ref. [13] in order to take non-adiabatic effects 
into account to correct the calculation of the tunneling probability in the adiabatic approximation (see 
ref. [ 12 ] for the derivation of these equations and the details.). 

Using the screening potential in the tunneling region thus obtained, we calculate the tunneling prob¬ 
ability in the presence of electrons by the following WKB formula 


P{E) = exp -2 


2M 


rRa 


I Rb 


dRJv(o){R) + AViR) - E 


exp y (iT[U(°) (i?) - E]^ exp J drAU(i?) 


( 11 ) 


where Ra and Rb are the classical turning points on both sides of the effective potential barrier (i?) -I- 
AV{R), and r^ and Tf, are the corresponding times along the imaginary time axis. We then convert the 
enhancement factor / = , where Pq{E') is the tunneling probability in the absence of electrons, 

into a screening energy using the relation 


= log f log 

nrj{E^) ^\Po{E')) 7Trj{E^) H Po{E]^) J 


( 12 ) 


where ri{E) is the Sommerfeld parameter, E^ is the kinetic energy of the relative motion between the 
target and projectile nuclei and e*-®! is the total energy of electrons in the center of mass system in the 
initial asymptotic region. The latter is identical with the screening potential AV at the initial time and 
is given by 


el®l = -f^eVT + er 


(13) 


. me~\-MT 
Mt 


where ex is the binding energy of electrons in the initial state in the target atom and vx = m^+Mp ’ 
Voo being aX t = —oo, is the velocity of the target nucleus relative to the center-of-mass of the 

projectile and target nuclei at the initial time. The reduced mass Ue is given by - 7 - = ^ 

f-ie ‘P'e 

Strictly speaking ex on the r.h.s. of eq.(13) should be replaced by e'j, = ex x ^ 

However, the difference between and ex can be ignored in the present problem. Note that we compare 
the tunneling probabilities for the same kinetic energy of relative motion of the nuclei in the presence 
and in the absence of the electrons. That is why we use different notations for the energy arguments in 
the barrier penetrability in the second term of eq.( 12 ). Also, we use in eq.( 12 ) and in what follows the 
lower index 0 to denote the barrier penetrability and the cross section calculated in a two body system 
and distinguish them from the corresponding quantities calculated including electrons. 
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We now apply our formalism to D+d and ^He“''+d reactions. We choose these systems for simplicity 
of the treatment because the screening effects are due to a single electron. Moreover, there exists 
experimental data for the D+d reaction [14] at a low energy, i.e. at Ec.m = 1-62 keV, though experiments 
have been performed for a molecular target rather than an atomic target [14]. 
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FIG. 1. Screening potential for the D+d reaction at the center-of-mass energies 1 and 200 keV as a function of 
the separation distance between the nuclei. The filled triangles show the position of the external classical turning 
point. 


Fig.l shows the screening potential for the D+d reaction at keV (solid line) and 200 keV 

(dashed line). The asymptotic values and their incident energy dependence can be understood from 
eq.(13). The closed triangles show the external classical turning points. Two interesting things can be 
noticed. The first is that the value of the screening potential at the external turning point for 1 keV is 
-34.0 eV, which matches with the average of binding energies e[fj^=-54.4 eV in the lowest gerade and 
e[](4=-13.6 eV in the ungerade states, i.e. in the Is- and 2p- states, of the united atom "‘He''". This 
indicates that the reaction takes place almost adiabatically in both gerade and ungerade configurations 
at this energy. The second observation is that the screening potential for Ec.m. = 1 keV changes 
very fast just inside the external classical turning point. This can be understood from eqs.(9) and (7) 
as a consequence that the contribution to the mean potential from the ungerade configuration, which 
has higher electronic energy, quickly dies out as the relative motion between the projectile and target 
penetrates into the tunneling region. 
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FIG. 2. The same as Fig.l, but for the ^He''’-|-d reaction at three incident energies Ec.m.= 2.5, 5.0 and 500 
keV. 


Fig.2 shows the screening potential for the ^He+ -|-d reaction. Though it is not so drastic as that 
for the case of D-l-d reaction, we see a similar structure near the external classical turning point for 
the cases of =2.5(the solid line) and 5.0 keV(the dotted line). The screening potentials for these 
two energies merge for distances smaller than the external classical turning point for Ec.m.=5.0 keV. In 
this system, the radial variation of the screening potential is caused by the admixture of excited states 
of electron due to pre-tunneling transitions, i.e. due to the transitions of the electronic state induced 
on the way from the initial asymptotic region to the classical turning point. One will then expect that 
the variation of the screening potential just after the system entered the tunneling region will get less 
significant at lower incident energies, because the process will become more and more adiabatic and the 
pre-tunneling electronic transition will become less important. This accords with Fig.2, where one sees 
that the variation of the screening potential at ifc.m.=2.5 keV is much smaller than that for ifc.m.=5.0 
keV. This contrasts to the case for D-l-d reaction, where the symmetry of the system admixes the gerade 
and ungerade states with equal weight at any incident energies including the low energy adiabatic limit. 
Note that eqs.(7) and (9) indicate that the screening potential at low energies will converge to the value 
in the adiabatic limit given by the binding energy of the Is state of electron in the united atom as the 
tunneling process proceeds, and thus explains why the screening potentials for ifc.m.=2.5(the solid line) 
and 5.0 keV(the dotted line) merge for distances smaller than the classical turning point for ifc.m.=5.0 
keV. Note that the screening potential at ifc.m.=500 keV (the dashed line) shows no corresponding sharp 
radial variation and looks quite different from that for ifc.m.=2.5(the solid line) and 5.0 keV (the dotted 
line). This will be partly because the tunneling process at =500 keV is near to the sudden limit 
rather than to the adiabatic limit, and also because the tunneling region is much smaller. Note also that 
the abscissa is log-scale. It makes the variation at larger values of R appear more drastic. The situation 
is similar in the case of D-(-d fusion reaction at ifc.m. =200 keV shown in Fig.l. 
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FIG. 3. Screening energy Ue as a function of the incident center-of-mass energy for the D+d reaction. The 
experimental value is for a molecular deuteron target taken from ref. 14. 
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Fig.3 shows the screening energy for the D+d reaction. The closed circles are the results of our method. 
The open squares have been calculated in the same way as in ref. [10]. The horizontal solid and dashed 
lines are the screening energies in the sudden and adiabatic reaction limits, and , respectively, 

which are given by 


[/f) = X 2 X ZpZt X en, 

® Mp + Mt 

= .. . X 2 X ZpZp X 13.6eV" = 13.6ey. 

Mp + Mt 

= i [(54.4 - 13.6) + (13.6 - 13.6)] 

= 20.4ey 


(14) 

(15) 

(16) 

(17) 

(18) 


These formulas, eqs.(14) and (16), for limiting cases are derived under the assumption that the tunneling 
region is much smaller than the Bohr radius of the united atom. Moreover, eq.(14) has assumed that the 
screening electron occupies the Is state of the target atom (see ref. [12] for details). In eq.(14), cr = 13.6 
eV is the binding energy of the Is orbit in the Hydrogen atom. In eq.(16), ep and cjja are the binding 
energies of the electron in the target and united atoms, respectively. The electron is assumed to occupy 
the adiabatic state with the same label i in both atoms because of the slow adiabatic process. In the 
second line of i.e. in eq.(17), we have used the actual values in the present case by taking the 

symmetry property of the D+d system into account. As one expects, the screening energy converges to 
that in the sudden reaction limit at high energies. It converges to the adiabatic limit at low energies if 
one calculates in the way of ref. [10] by studying only the classical region. The star with error bar is the 
experimental value taken from ref. [14]. However, this should be taken merely as a reference, because 
as mentioned before the experiments have been performed not for an atomic deuteron target, but for 
a molecular deuteron target. In this connection, we wish to mention that ref. [11] has shown that the 
screening effect for the molecular D 2 target is larger than that for the atomic D target and attributed 
this difference to the reflection symmetry of the d+D system. 

The remarkable thing is that our calculations give systematically a larger screening energy than that 
in the conventional calculations. At low energies, this can be understood in the following way. Using 
the screening potential at the external turning point Rt, the enhancement factor is calculated in the 
conventional method, e.g. in ref. [10], by 

, ao(U^+e«-AU(A-t)) 

A - -^ 
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( 20 ) 


_ CJo{E^ +€T-^V{Rt)) 

MEr) 

_ ao(£;g + + £/i"V2) 

In transforming from eq.(19) to eq.(20), we have ignored the difference between and ex given by 
eq.(13) in accord with the adiabatic process. Also, in order to move further to eq.(21), we have used the 
fact, which we remarked before concerning Fig.l, that the screening potential at the external classical 
turning point can be understood in terms of the binding energies of the electron in the gerade and 
ungerade configurations of the united atom. On the other hand, our method, which handles the tunneling 
region explicitly, leads to 


_ cr(9) + cr(“) 

~ ^0 

^ ao(Ag + + ao(Ag + 

2ao{E^) 


( 22 ) 

(23) 


for the enhancement factor. These equations can be derived from eqs.(7),(9) and (11) by assuming that 
there are no change with the adiabatic energies and the adiabatic states in the tunneling region. This 
assumption will be reasonable in the present case, for which the tunneling region is much smaller than 
the Bohr radius of the united atom as shown in Table 1. Since the excitation function of the fusion cross 
section is a convex increasing function of the incident energy, ft is larger than fc- The conventional 
method thus underestimates the screening energy. 

Fig.3 clearly exemplifies this effect. It is important to properly calculate the enhancement factor in 
order to get a reliable value of the screening energy. This can be achieved either by explicitly handling 
the tunneling region like in our method, or by studying the distribution of the electronic state over 
different adiabatic states at the external classical turning point, and calculate the fusion probability for 
each of them and taking average afterwards with the proper weight. In this respect, we note that our 
numerical result at iJc.m.=1.6 keV agrees with the value given by eq.(23). Assume that the electronic 
wave function is distributed over different adiabatic states (/?„, whose corresponding binding energy 
et — Un, with probability P„, then eqs.(7), (9) and (11) lead in general to 


„ _ 'EnPn X ^o{E^ + Un) 

ME]!) 

which should be compared with the formula 

r _ ^o(Pk a EnPn X Un) 

MEf) 


(24) 


(25) 


in the conventional method, where one first calculates the average potential and then calculates the 
tunneling probability for it. We have assumed that the tunneling region is much smaller than the Bohr 
radius of the united atom and ignored the change of the adiabatic states and their corresponding energies 
in the tunneling region. Note that the presence of the denominator on the r.h.s. of eq.(9) is essential in 
deriving eq.(24). 


TABLE 1. Incident energy dependence of the external classical turning point for the D+d reaction. The initial 
velocity of the relative motion in the unit of the light speed c, Voo/c, the external turning point in the absence 
of an electron and in its presence Rt, and the mean square radius of the electron at the external classical 
turning point are shown for three bombarding energies. 


Px(keV) 

Voojc 

A™ (A) 

At (A) 

Electron Radius(A) 

100.0 

0.0146 

0.000144 

0.000144 

0.456 

10.0 

0.00462 

0.00144 

0.00144 

0.518 
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1.0 


0.00146 


0.0144 


0.0142 


0.688 


In summary, we have presented a semiclassical mean field theory of quantum tunneling which treats 
both classical and tunneling regions in a consistent way. Applying the formalism to the problem of 
screening effects by bound target electrons in low energy nuclear reactions in laboratories, we have 
shown that the screening potential shows a characteristic radial variation contrary to the assumption 
of a constant potential shift in all previous analyses. We have shown also that the proper treatment of 
the tunneling region leads to an increase of the screening energy compared with that estimated in the 
previous mean field theory, which studies only the classical region and calculates the tunneling probability 
by using the average potential at the external classical turning point. The above effects are, however, 
too small to explain the large experimental screening energies reported in ref. [2]. Remember in this 
connection that the large screening energy obtained in this study for the D+d reaction has been caused 
by the symmetry special to this system, and cannot be generalized to other systems. 
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